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Abstract 

In this article, we study the radiative transitions among the vector and scalar heavy quarko- 
nium states with the covariant light-front quark model. In calculations, we observe that the 
■ radiative decay widths are sensitive to the constituent quark masses and the shape parameters 

of the wave-functions, and reproduce the experimental data with suitable parameters. 
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• Recently, the BESIII collaboration observed the hrst evidence for direct two-photon transition 

II | ip' — > J/ipj^f with the branching fraction (3.3 ± 0.6j^f) x 10 -4 in a sample of 106 million ip' 

Qh decays collected by the BESIII detector [T|. The branching fractions of the double E\ transitions 

ip' — > (J/ipj) Xc jJ through the intermediate states Xcj (j = 0,1,2) are also reported Q], while 
previous experimental data indicates that the double radiative decays ip' — > J/ijj-f-f take place 
through the decay cascades ip' —> Xcjl ~~ ^ J I^Pll with tiny non-resonance contributions [H In 
Ref.[4], He et al study the discrete contributions to decays ip' — > J/ip~fy due to the E\ transitions 
using the heavy quarkonium effective Lagrangian [5]. No theoretical work on the non-resonance's 
qq ■ contributions exists. On the other hand, we expect that there are non-resonance's contributions to 

\Q ' the doubly radiative decays among the bottomnium states T(nS), the doubly radiative transition 

■ T" — > T'77 has been observed [5]. 

\q j The radiative transitions among the heavy quarkonium states are usually calculated model- 

dependently by the nonrelativistic potential quark models with considerable relativistic corrections 
[7] , or calculated model-independently by the lattice QCD [8] and effective field theory [5] , one can 
consult the recent review |10] for more references. In general, we expect to study both the reso- 
nance's and non-resonance's contributions in the doubly radiative transitions with the covariant 
light-front quark model, where the wave-functions are expressed in terms of the internal variables 
of the quark and gluon, and maintain Lorentz covariance. In Refill], Jaus introduces the covari- 
ant light-front quark model and takes into account the zero-mode contributions systematically to 
preserve covariance and remove dependence of physical quantities on the light-front direction. The 
covariant light-front quark model has been successfully applied to calculate the S-wave and P-wave 
meson's decay constants and form- factors [HI [T3J 03] ■ In Refs. |151 116) . the authors study the Mi 
transitions T(nS) — > ^(n'S) 7 in the covariant light-front quark model, and observe that the Mi 
transitions are sensitive to the heavy quark masses and shape parameters of the wave- functions, 
the existing experimental data cannot be reproduced consistently with suitable parameters |15j . 
It is interesting to study whether or not the covariant light-front quark model can be successfully 
applied to calculate the E\ transitions among the heavy quarkonium states. 

In this article, we intend to study the radiative transitions (the E\ transitions) among the 
vector and scalar heavy quarkonium states using the covariant light-front quark model as the 
first step, because they can be taken as the sub-processes of the doubly radiative decays, and the 
calculations are relatively simple. Experimentally, the widths of the radiative decays ip' — ¥ Xcoli 
XcO — > J/^Pl: T' — > Xbol, T" —> Xbol, T" — > Xwl nave been precisely measured [6]. A large 
amount of bottomonium states will be produced at the Large Hadron Collider, and the radiative 
transitions will be studied experimentally. We can study the radiative transitions by measuring 
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Figure 1: The Feynman diagrams contribute to the form- factors. The photon is emitted from the 
quark (antiquark) line in the diagram A (B). 



the energy spectrum of the photons or reconstructing the final quarkonium states, although the 
soft photons are difficult to identify. The four-photon decay cascades T(35) —> -fXbj(2P) 
77T(1D) — > 777X6j(l-P) —> -f-fjjT(lS) —> 7777^+^" (j = 0,1,2) have been observed by the 
CLEO collaboration [T7], where the softest photons have the energy less than 90MeV. 

The article is arranged as follows: we calculate the radiative transitions among the vector and 
scalar heavy quarkonium states with the covariant light-front quark model in Sect. 2; in Sect. 3, we 
present the numerical results and discussions; and Sect. 4 is reserved for our conclusions. 



2 Radiative decays with covariant light-front quark model 

The radiative transitions among the heavy quarkonium states can be described by the following 
electromagnetic lagrangian C, 



L = 



-eebbj^bA* 1 - eecCj^cA* 1 . 



(1) 



where the A^ is the electromagnetic field, the e is the electromagnetic coupling constant, 



and 



The transition amplitudes £ M 7~J /_i '' s ' 7 and £ fl A^ v ' y can be decomposed as 



V 



fv Sl (e^q ■ P> - P;e ■ q) , 

gsvj (e^9 • p - P^* ■ q) , 



(2) 



according to Lorentz covariance, where the V and S denote the vector and scalar mesons respec- 
tively, the P, P' , q are the momenta of the initial mesons, final mesons and photons respectively, 
the £ M and are the polarization vectors of the photons and vector mesons respectively, the fvs-y 
and gsv-y are the electromagnetic form-factors (or the coupling constants) at q 2 = 0. The am- 
plitudes T^^ Sl and A s ^ Vl satisfy conservation of electromagnetic currents, we can replace the 
photon's polarization vector £ M with its momentum q^ in the transition amplitudes ^T^~ 7 and 
f y A S-*V 1 tQ btain q^r^ Sl = q^A^ v i = 0. 

From the lagrangian C, we can draw the Feynman diagrams (see Fig.l) and write down the 
transition amplitudes, 
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the Hy and Hg are the light-front wave-functions of the vector and scalar mesons respectively, the 
A and B denote the diagrams in which the photon emitted from the quark and antiquark lines, 
respectively (see Fig.l), and the Wy is a parameter. In writing the transition amplitudes T^^ Sl 
and A^ Vl , we have used the following definitions of the quark- meson-antiquark vertexes iYm 



iT 



M 



-iHc 



7m 



(Pi - P2) A 



W v 



for scalar meson , 
for vector meson . 



(7) 



where the p\ and P2 are momenta of the quark and antiquark, respectively. 

We assume that the covariant light-front wave-functions Hy and Hs are analytic in the upper 
(lower) complex plane, close the integral contour in the upper (lower) p~[ complex plane for the 
diagram A (£?), which corresponds to set the antiquark (quark) on the mass-shell. The one-shell 
restrictions are implemented by the following replacements, 



P2 

Ni 
Hy 



P? 



d A pi 



for the Feynman diagram A and 



-> P2 = m Q I Pi -> Pi = P - P2 

-> #i=p?-m£, N[^N[ = 
-» Hy = hy(x2,p±) , Hs -> Hs = h s {x 2 ,px) 

dx 2 d 2 p± 

x 2 N 1 N[ ' 



Wy — > Wv = Wy , 



(8) 
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m Q, P2 -> f>2 = P - Pi , 

-> N 2 =p 2 2 -m 2 Q , N£-> N^ = pf - m 2 Q , 

-» H v = h v (x 2 ,p±) , H s -> H s = h s (x 2 ,px) 
dx 2 d 2 p± 
x^N' ' 
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(9) 



for the Feynman diagram S, where we have used the light-front decomposition of the momenta 

P = (P+,P-,0 X ), P ± = P°±P\ pf 2 = X 1)2 P + , Xl +X 2 = 1, Pl,2X = ±PX, 9 = (0,?-,«±) = 
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(0, q , 0j_). Then we obtain the form-factors or coupling constants at zero momentum transition^), 
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where the light-front wave-functions are defined as 
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2 For technical details, one can consult Refs.[TTlll2|. 
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the Mi and Mf are the masses of the initial and final heavy quarkonium states, the can be 
viewed as the energy of the quark (antiquark), and the Mq can be viewed as the kinetic invariant 
mass of the quark- antiquark system. 

In calculations, we have used the following rules [TT| IT2] . 

pj* = P^A^ + if , 

ft ft = g^Af* + P' 1 P v Af > + (P"g" + q^P u )A (2) + q^q v A {2) , 

N 2 -> Z 2 , 

/,V^ - -i/^ P i,/^^(^-2^>)=0, 

_ X l 4 (1) _ Pl'gl 4 (2) _ 2 (P± ■ g-L) 2 

4 (2) 

4(2) _ 4(1) Al) 4 (2) _ .(1) ,(1) ,(2) _ A1) A (l)_ £l_ n „x 
— ^ A] , — /i x ^i 2 > Aj ~ ^2 ^2 ^2 ' \ ' 

In the covariant light-front quark models, we usually take the three-dimensional harmonic- 
oscillator wave-functions in momentum space to approximate the radial wave-functions [121 1151 
ITS] . it is a phenomenological approach. The heavy quarkonium spectra are usually obtained by 
numerically solving the Schrodinger equation with the Cornell potential or other phenomenological 
potentials, such as the logarithmic potential, the power-law potential, the QCD-motivated potential 
[19] . It is difficult to translate the numerical wave- functions into the momentum space and re- 
formulate them with the light-front variables. While the analytical three-dimensional harmonic- 
oscillator wave-functions are simple, and it is easy to translate them into the momentum space. 
Furthermore, we cannot solve the light-front hamiltonian as an eigenvalue problem analytically to 
obtain the light-front wave-functions. 



3 Numerical results and discussions 

We take the masses of the heavy quarkonium states from the Particle Data Group, M Xe0 = 
3414.75 MeV, Mj/j, = 3096.916 MeV, M$, = 3686.09 MeV, M T = 9.46030 GeV, M r > = 10.02326 GeV, 
M T " = 10.3552 GeV, M Xb0 = 9.85944 GeV, M x > w = 10.2325 GeV [6], and obtain the radiative de- 
cay widths 

r Xe0 ->J/*7 = 2214.63 e 2 c e 2 g 2 XcoJ/qp7 KeV , 
I>_> Xe07 = 473.535 e 2 e 2 f 2 , Xc0l KeV , 
rr'-* M7 = 113.783 e 2 b e 2 f r , Xbol KeV , 
r T ^ XM)7 = 3005.49 e 2 e 2 f 2 „ Xb0l KeV , 

IV^ o7 = 48.135 e 2 e 2 f 2 „ x , o7 KeV , (14) 
Compared to the experimental data from the Particle Data Group [B] , 

T Xc0 ->J/in = 122.85 ± 9.36 ± 8.4 (121.68 ± 7.02 ± 8.32) KeV , 

rV_>. Xo07 = 27.6848 ± 1.5488 ± 0.8866 (29.4272 ± 0.8712 ± 0.9424) KeV , 

r T'^ Xb07 = 1.21524 ±0.09994 ±0.12792 KeV, 

r x „^ Xbo7 = 0.06096 ±0.004995 ±0.008128 KeV, 

r T"-> x ' 60 7 = 1-19888 ±0.10915 ±0.12192 KeV, (15) 

where the first uncertainties come from the total widths and the second uncertainties come from 
the branching ratios, we can draw the conclusion tentatively that the coupling constants fyg~ 
differ greatly from each other for the bottomonium states, while g 2 x ji^ 1 ~ f%>'xeoT 
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From Eqs. (10-12), we can see that the form-factors or coupling constants depend on two kinds 
of inputs parameters, the constituent quark masses ttiq and the shape parameters /3 of the wave- 
functions. The spin averaged ground state masses are m "c+3 m j/v _ 3. 068 GeV, — — + ^ Mt — 
9.443 GeV from the Particle Data Group [6], we estimate the constituent quark masses m c ~ 
1.5 GeV and nib ~ 4.7 GeV. In this article, we take the constituent quark masses tuq and shape 
parameters (3 of the wave-functions as free parameters and search for the optimal values at the 
ranges m c = (1.3 - 1.7) GeV and m b = (4.5 - 5.0) GeV. 

In numerical calculations, we observe that the form-factors (or coupling constants) are sensitive 
to the constituent quark masses tuq and the shape parameters j3, small variations of those input 
parameters can lead to large changes of the form- factors (thereafter the decay widths) ; the optimal 
values are m c = 1.5 GeV and m b = 4.8 GeV, (3 Xc0 = 1.30 GeV, /3 J/V , = 0.76 GeV, fa, = 0.83 GeV, 
P r , = 1.30 GeV, p r » = 0.976 GeV, /3 Xb0 = 0.94 GeV, /3 x - w = 0.945 GeV, which can reproduce the 
decay widths of the five observed processes \c0 — > J/tpJ, "0' ~~ ► Xco7j T' — > XM)7, T" — > Xf>o7 and 
Y" Xf,o7 0- ^ the form- factors (or coupling constants) are not sensitive to the shape parameters 
/3 of the wave-functions, we can introduce two parameters to characterize the charmonium and 
bottomonium states respectively, and search for the optimal values to reproduce the decay widths 
of the five observed radiative transitions, however, two parameters cannot lead to satisfactory 
results. 

The optimal values m c = 1.5 GeV and mj = 4.8 GeV approximate or equal to the estimated 
values m c rj 1.5 GeV and m& ~ 4.7 GeV. In the constituent quark model, the usually used 
constituent quark masses are m c = (1.4 — 1.6) GeV and to& = (4-7 — 4.9) GeV, it is natural to take 
the values m c = 1.5 ± 0.1 GeV and m b = 4.8 ± 0.1 GeV. 

We choose (or suppose) the shape parameters have the same uncertainties as the constituent 
quark masses tentatively, i.e. /3 Xc0 = 1.30 ± 0.1 GeV, /3 J/4 , = 0.76 ±0.1 GeV, /V = 0.83 ±0.1 GeV, 
Pr> = 1.30 ± 0.1 GeV, /3 T " = 0.976 ± 0.1 GeV, /3 Xb0 = 0.94 ± 0.1 GeV, /3 Xbo = 0.945 ± 0.1 GeV, and 
study the resulting uncertainties, 

r — 191 r:/i+160.22 +34.32 +39.89 Tf„v 

1 XcO^J/4>7 ~ lzl - d ^-97.18 -26.26 -32.67 - tYL v ' 

F — 07 qh+26.67 +9.83 +1.56 tv- Tr 

1 0'^Xco7 — z '- yu -17.97 -9.10 -3.55 lvev ' 

p _ i 2 o+2.31 +0.07 +0.52 K V 

1 T'->Xbo7 — 13 -0.13 -0.34 ■ rve v > 

r~ n nfinR+ 0078 +0.2160 +0.2343 K v 

L ^"^XbOl ~ u - uouo -0. 0072 -0.0605 -0.0589 v ' 

rx-^ o7 = i.i9±g:g±;:fg±?:SKev, (ie) 

where the uncertainties come from the heavy quark masses and shape parameters of the initial and 
final quarkonium states, sequentially. 

From Eqs. (15-16), we can see that the central values of the present predictions are consistent 
with the experimental data from the Particle Data Group [6] , however, the uncertainties come from 
the heavy quark masses exceed 100% for the transitions Xco - > J/^li ip' ™>> Xco7, T' — > Xbo7, and 
the uncertainties come from the shape parameters also exceed 100% for the transitions T" — > Xbol, 
T" — > x'bol- We have little room for varying those parameters to reproduce the experimental data, 
which weakens the predictive ability of the covariant light-front quark model remarkably. On the 
other hand, if we take the uncertainties of the experimental data on the radiative transitions (see 
Eq.(15)) to constrain the uncertainties of the constituent quark masses and shape parameters, the 
allowed uncertainties are \Sm c \ < lOMeV, \Sm b \ < lOMeV, |<5^ Xc0 | < 35MeV, \5/3 J/4 ,\ < 30MeV, 
|(5/3. '| < 17MeV, \6/3 r >\ < 100 MeV, \6/3 r »\ < 8MeV, \S/3 Xw \ < lOMeV, \&P X >J < 6MeV. 

The heavy quarkonium states have equal constituent quark masses, niQ — m-Q, the heavy quark 
masses are taken as one parameter niQ rather than two parameters mg and m-Q, see Eqs. (10-12), 
the total uncertainties come from the heavy quark masses are larger than that come from the heavy 
quark masses ttiq and m-Q, respectively. 

In Figs. 2-3, we plot the radiative decay widths r Xc0 ^j/^, 7 , r^,/_ >Xc07 , Tr'~^ Xwn , Tr"^ Xb0 -y and 
rY"->x 60 T with variations of the shape parameters /3. From the figures, we can see that the widths 
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are sensitive to the shape parameters /3 indeed, small variations of the shape parameters /3 can 
lead to rather large changes of the decay widths, especially when there are nodi in the radial 
wave-functions, it is very difficult (or impossible) to obtain two universal parameters Pec and /3^ b 
to characterize the charmonium and bottomonium states respectively. 

We can take the optimal values m c — 1.5 GeV and nib = 4.8 GeV, (3 Xc0 = 1.30 GeV, /Jjm = 
0.76 GeV, /V = 0.83 GeV, /3 T < = 1.30 GeV, /3 T » = 0.976 GeV, /3 XM = 0.94 GeV, /3 x > o = 0.945 GeV 
as the basic input parameters and study other radiative transitions among the heavy quarkonium 
states. For example, we can take the parameter f3 Xb0 — 0.94 GeV to calculate the decay width 
r XM) ^T7- From Fig. 3, we can see that the prediction r Xb0 ^x7 = (18 — 34) KeV for (3y = (0.4 — 
1.4) GeV is consistent with the value r Xb0 _ > x7 = (24—26) KeV from the potential quark models [2D]. 
If we take the parameter /3x = 1.16 GeV, the predictions T Xi0 ^~cy — 24.75 KeV and T x > ^.y^ = 
5.97 KeV are consistent with the values from the potential quark models r Xb0 ^x7 = (24 — 26) KeV 
and r x ' bo _i.x7 = (4.5 — 8.5) KeV [20], while the prediction r Xtf ^x'7 = 46.74 KeV is much larger 
than the value r^^x'7 = (H — 12) KeV from the potential quark model [20] . 

Without precise experimental data, we cannot determine the value of the shape parameter /3x- 
On the other hand, there are controversies for the spectroscopy of the charmonium states, and lack 
experimental data on other radiative transitions among the vector and scalar charmonium states. 
We prefer study those processes in the future. 

In Ref.[15j. W. Wang calculates the S'-wave to S'-wave radiative transitions T — > r/bj, T' — » 
rj b -f, V -> 77^7, T" -> r\bl, T" -> 77^,7, J/tp -> r] c j, ip' -> ri c j, ip' -> 77^,7 with the covariant 
light-front quark model, and observes that those Mi transitions are sensitive to the heavy quark 
masses and shape parameters of the wave-functions, the existing experimental data cannot be 
reproduced consistently with suitable parameters. The input parameters used in studying the 
observed Mi transitions again fail to reproduce the experimental data on the Ei transitions of 
the heavy quarkonium states. In Ref. 21 , H. W. Ke et al observe that the three-dimensional 
harmonic-oscillator wave-functions lead to obvious contradictions between the experimental and 
theoretical values of the T(nS) decay constants. They modify the wave- functions by introducing 
several additional parameters to describe the radial quantum number n dependence and keep the 
orthogonality among the nS states, then fit the decay constants to determine the T(nS) wave- 
functions and study the radiative decays T(nS) — > rjj, + 7, the prediction of the T(3S f ) — > rjb + 7 
is in accordance with the experimental data by the order of magnitude, the discrepancy between 
the experimental and theoretical values is not well smeared. In the present case, the different 
values of the shape parameters /3 Xc0 , /3j/^>, /3^>, /3x', Px", Pxboi Px' b0 violate the orthogonality of 
the wave-functions, it is the shortcoming of the present work. We can also draw the conclusion 
tentatively that we cannot get rid of the sensitivity to the shape parameters /3 without introducing 
additional parameters in the present case, however, the experimental data are far from enough to 
fit the additional parameters. 

4 Conclusions 

In this article, we study the radiative transitions among the vector and scalar heavy quarkonium 
states in the framework of the covariant light-front quark model. In calculations, we observe that 
the radiative decay widths are sensitive to the constituent quark masses and the shape parameters 
of the light-front wave-functions. We reproduce the experimental data for the observed processes 
with suitable parameters, while the predictions for the un-observed processes are consistent or 
inconsistent with other theoretical calculations. The parameters can be fitted to the precise ex- 
perimental data in the future, although it is a difficult work. 
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Figure 2: The radiative decay widths with variations of the shape parameter /3 Xc0 . In (I) Xco 
J/ipj, A, B, C, D correspond to pj/^/P Xc0 = 0.56/1.3, 0.66/1.3, 0.76/1.3, 0.86/1.3, respectively; 
in (II) ip' -> Xcol, A, B, C, D correspond to fy>//3 Xc0 = 0.63/1.3, 0.73/1.3, 0.83/1.3, 0.93/1.3, 
respectively. 




Figure 3: The radiative decay widths with variations of the shape parameters /3 Xb0 or /3 X ' ■ In (I) 
T' Xbo7, A, B, C, D correspond to /3r'/P Xb0 = 1.1/0.94, 1.2/0.94, 1.3/0.94, 1.4/0.947 respec- 
tively; in (II) T" -> X607, A, B, C, D correspond to /3 T » / 'P X w = 0.776/0.94, 0.876/0.94, 0.976/0.94, 
1.076/0.94, respectively; in (III) T" ->■ x'tol, A , B , C, D correspond to /3 T » / ' P x > b0 = 0.776/0.945, 
0.876/0.945, 0.976/0.945, 1.076/0.945, respectively; in (IV) x&o T 7 , j3y = 0.94 b °GeV. 
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